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FOREWARD 

The  Investigation  describee  In  this  First  Yearly  Report  to  the 
Office  of  Naval  Research  (ONR)  Is  the  outcome  of  a  continuing  collabora- 
tlve  effort  between  the  respective  principal  Investigators,  Professor 
Joseph  A.C.  Humphrey  at  the  University  of  California,  Berkeley  (UCB), 
and  Professor  Brian  E.  Launder  at  the  University  of  Manchester  Insti¬ 
tute  of  Science  and  Technology  (UMIST).  Because  the  first  six  months 
of  this  work  were  conducted  while  Professor  Launder  was  at  the  University 
of  California,  Davis,  and  because  the  first  year  of  research  has  Involved 
an  especially  close  collaboration  between  the  principal  Investigators 
In  order  to  ensure  continuing  complementary  research  efforts  at  their 
respective  institutions,  it  Is  appropriate  that  the  first  yearly  report 
should  be  co-authored. 

At  UCB  an  experimental  facility  has  been  constructed  which  allows 
the  measurement  of  fluid  velocity  and  turbulence  characteristics  In  and 
downstream  of  curved  ducts  of  square  or  circular  cross-section  using  an 
automatic  data  acquisition  system.  A  corresponding  heat-transfer  experi¬ 
ment  for  flow  in  ducts  of  circular  cross-section  is  In  preparation  at  UMIST. 
A  numerical  procedure  for  predicting  turbulent  flow  In  and  downstream  of 
curved  ducts  of  square  cross-section  is  under  development  at  UCB.  A 
corresponding  effort  for  predicting  flows  in  ducts  of  circular  cross- 
section  is  being  developed  at  UMIST.  It  Is  Important  to  note  that  while 
the  two  calculation  schemes  differ  in  their  coordinate  geometry  (cylindri¬ 
cal  versus  toroidal)  the  numerical  approach  and,  especially,  the  turbulence 
modeling,  are  the  same.  As  work  progresses,  measurements,  calculations  and 
the  numerical  schemes  under  development  at  the  two  Institutions  will  be 
critically  compared  and  evaluated  for  the  benefit  of  both. 
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In  several  Important  ways  the  first  year  of  research  for  ONR  has 
overlapped  with  a  one  year  program  of  research  conducted  by  Professor 
Humphrey  at  UCB  for  the  National  Science  Foundation  (NSF).  Thus,  while 
the  ONR  program  has  benefited  from  the  experimental  facility  and  pre¬ 
liminary  numerical  explorations  Initiated  with  NSF  funds,  the  NSF  pro¬ 
gram  has  benefited  from  continuing  and  extensive  numerical  developments 
made  possible  by  ONR.  We  believe  that  this  overlapping  of  programs, 
from  which  a  mere  duplication  of  research  efforts  has  been  strlckly 
barred,  has  strongly  benefited  both  Investigations  and  the  study  of 
curved  duct  flows  In  general.  We  wish  to  express  our  appreciation  to 
Dr.  Win  Aung  of  the  National  Science  Foundation  and  Dr.  Keith  Elllngsworth 
of  the  Office  of  Naval  Research  for  permitting  and  encouraging  this 
mutually  beneficial  program  interaction. 

This  report  is,  primarily,  an  account  of  the  work  performed  for  ONR 
throughout  the  first  year  (1979-1980)  of  research.  However,  because  the 
work  overlaps  with  the  study  performed  for  NSF  during  the  same  period  of 
research,  reference  will  be  made  to  the  latter  and  to  an  oncoming  final 
NSF  report.  Upon  completion,  the  NSF  document  will  serve  as  an  iterlm 
ONR  report  and  will  be  made  available  to  those  Individuals  on  the  ONR 
distribution  list  wishing  to  receive  a  copy. 
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ABSTRACT 

An  account  Is  rendered  of  the  research  conducted  during  the  first 
year  of  an  Investigation  aimed  at  developing  a  computational  facility 
for  predicting  turbulent  flow  in  180°  curved  passages  of  rectangular 
and  circular  cross-section.  In  order  to  evaluate  critically  the 
numerical  scheme  and  turbulence  model  under  development  an  experimental 
apparatus  for  measuring  the  velocity  and  turbulence  characteristics  of 
such  flows  has  been  constructed. 

The  first  year  of  research  has  been  dedicated  primarily  to: 

1.  Developing,  testing  and  running  the  experimental  apparatus  and 
associated  electronic  Instrumentation  (Including  a  computer- 
aided  automatic  data  acquisition  system). 

2.  Developing  and  testing  a  quadratic-upstream  convective  dif¬ 
ferencing  scheme  to  alleviate  the  problems  of  numerical  or 
"false"  diffusion  associated  with  simple  upwind  differencing. 

Progress  made  on  these  two  accounts  is  reviewed  particularly  in  relation 
to  the  overall  objectives  of  this  research  program.  An  evaluation  Is 
made  of  the  present  status  of  the  program  and  a  tentative  schedule  of 
projected  activities  is  outlined  for  years  two  and  three  of  research. 


1.  REMINDER  OF  THE  PROBLEM  OF  INTEREST 


It  has  only  been  relatively  recently  (the  past  6-7  years)  that 
numerical  techniques  and  models  of  turbulence  have  been  applied  to  the 
problem  of  predicting  complex  three-dimensional  (3-D)  flows  with  heat 
transfer.  While  there  is  a  clear  indication  in  the  literature  that 
some  measure  of  success  has  been  achieved  [1-10]  there  remain  serious 
uncertainties  relating  to: 

a)  The  efficiency  and  accuracy  of  the  various  numerical  pre¬ 
diction  schemes. 

b)  The  accuracy  and  generality  of  application  of  the  various  model 
of  turbulence. 

c)  The  rlgurousness  and  generality  of  the  tests  applied  to  the 
numerical  schemes  and  models  of  turbulence  respectively. 

Point  (c)  above  is  especially  important  since  the  rigorous  testing 
of  a  turbulence  model  requires  that  numerical  aspects  of  the  problem  be 
error  free.  Thus,  for  example,  it  would  be  desirable  to  reduce  to  a 
minimum  numerical  or  "false"  diffusion  in  a  numerical  procedure  prone  to 
this  effect  prior  to  establishing  the  relative  merits  (or  demerits)  of  a 
particular  model  of  turbulence. 

The  present  work  addresses  each  of  the  above  three  points  by  refer¬ 
ence  to  3-D  turbulent  flows  within,  and  downstream  of,  180°  curved  pas¬ 
sages  of  circular  and  square  cross-section.  The  choice  of  flow  was 
dictated  by  two  major  considerations: 

1.  That  is  arise  In  the  components  of  major  industrial  processes, 
such  as,  for  example,  bends  In  heat  exchangers,  and  hence  be 
of  practical  Importance. 
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2.  That  Is  permit  a  rigorous  and  fairly  general  test  of  models 
of  turbulence  purporting  to  predict  complex  3-D  flows. 

The  flow  entering  a  curved  passage  experiences  a  radial  pressure 
gradient-centrifugal  force  imbalance  which  subjects  slow-moving  fluid 
near  the  side  walls  (parallel  to  the  plane  of  curvature)  to  a  pressure 
force  which  drives  the  flow  from  the  outer  radius  wall  towards  the  inner 
radius  wall.  Continuity  demands  that  a  corresponding  flow  be  set  up 
along  the  curved  passage  symmetry  plane  from  the  inner  to  outer-radius 
wall.  This  motion  is  referred  to  as  secondary  motion  of  the  "first"  kind. 

In  the  straight  section  following  the  curved  passage  the  secondary  motion 
relaxes  through  turbulence  diffusion  and  redistribution  processes.  If 
the  straight  duct  is  long  enough  this  cross-stream  flow  is  replaced  by  a 
weaker  turbulent  stress-driven  secondary  motion  of  the  "second"  kind. 

The  transition  between  the  two  states  of  cross-stream  motion  can 
only  be  predicted  by  a  turbulence  model  in  which  Reynolds  stresses  are 
calculated  from  their  respective  transport  equations.  In  addition,  the 
model  must  be  able  to  account  for  outer-wall  destabilizing  effects  and 
inner-wall  stabilizing  effects  on  the  flow  in  die  curved  section.  It  is 
clear  that  the  combination  of  a  curved  duct  followed  by  a  straight  section 
represents  a  severe  and  very  general  test  case  for  any  numerical  model 
of  turbulence. 


2.  RESTATEMENT  OF  RESEARCH  PROGRAM  OBJECTIVES 


The  research  objectives  of  this  investigation  are  two: 

1.  To  develop,  test  and  apply  a  numerical  calculation  procedure 
capable  of  modeling  fairly  general  and  complex  3-D  flows,  strongly 
affected  by  turbulence  diffusion  and  redistribution  processes. 

2.  To  obtain  the  necessary  experimental  information  by  reference  to 
which  a  precise  and  fairly  extensive  evaluation  of  the  numerical 
model  may  be  carried  out. 

With  respect  to  the  first  objective,  attention  has  been  focused  on 
the  development  of  a  so-called  "semi -elliptic"  type  of  calculation  scheme. 
The  scheme  is  of  the  global  iteration  class  [11]  wherein  pressure  is  the 
only  primitive  variable  stored  three-dimensionally.  A  forward  marching 
algorithm  is  used  as  part  of  a  global  iteration  process  in  which  stream- 
wise  diffusion  contributions  to  the  balance  equations  are  neglected. 
Although,  due  to  the  nature  of  the  forward-marching  approach,  imbalances 
(in  continuity  or  in  the  representation  of  pressure  gradients)  are  intro¬ 
duced  in  the  system  of  governing  equations,  the  iterative  nature  of  the 
procedure  allows  their  reduction  by  systematic  correction  until  they  are 
reduced  to  a  negligible  level. 

The  present  scheme  is  based  on  a  3-D  elliptic  numerical  procedure 
developed  oy  Humphrey  [8,12]  but  is  being  modified  here  to  include  the 
semi-elliptic  method  of  Pratap  [13]  and  the  QUICK  scheme  for  convective 
differencing  of  Leonard  [14], 
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The  second  objective  has  been  dictated  by  the  general  unavailability 
of  sufficiently  extensive  and  accurate  experimental  data  concerning  the 
turbulence  characteristics  of  complex  3-D  flows.  The  laser-Doppler 
velocimeter  data  measured  by  Humphrey  [12]  (see  also  reference  [15]) 
alleviates  the  situation  somewhat  in  that  is  provides  a  test  for  the 
numerical  prediction  of  turbulent  flow  development  in  90°  curved  ducts. 

However,  to  our  knowledge  there  is  no  equivalent  body  of  information  for 
a  180°  curved  passage  connected  to  a  straight  downstream  section,  in 
which  the  pressure  dominated  flow  leaving  the  bend  section  relaxes 
through  vigorous  turbulent  interactions.  An  experimental  apparatus  meeting 
the  needs  for  obtaining  such  information  with  the  help  of  a  PDP  11/34 
mini-computer  is  described  in  an  Appendix  attached  to  this  report. 
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3.  SUMMARY  OF  ACCOMPLISHMENTS  FOR  YEAR  ONE 

During  the  first  year  research  efforts  have  been  concentrated  on: 

1.  Building  and  testing  the  experimental  facility  for  measuring  develop¬ 
ing  turbulent  flows  in  ducts  of  circular  and  square  cross-section. 

2.  Developing  and  testing  a  quadratic-upstream  interpolation  scheme 
for  the  finite  differencing  of  convection  terms  in  the  Navier- 
Stokes  equations. 

The  experimental  facility  is  described  in  detail  in  the  Appendix. 

Since  its  completion,  however,  it  has  been  modified  in  several  respects. 

The  most  significant  changes  have  consisted  in: 

1.  Replacing  4  out  of  5  perforated  plastic  plates  in  the  upstream 
duct  by  3  metal  screens  of  graduated  mesh  sizes. 

2.  Removing  the  globe  valve  in  the  upstream  duct. 

3.  Removing  the  venturi  meter  from  the  downstream  duct  location  and 
placing  it  on  the  pump  line  feeding  the  constant  head  tank. 

In  combination,  these  modifications  have  led  to  a  much  more  uniform 
flow,  of  reproducible  turbulence  characteristics,  in  the  straight  section 
proceeding  the  180°  curved  section. 

The  experimental  system  has  been  tested  carefully  and  extensively.  At 
the  time  of  writing  measurements  are  being  completed  for  two  velocity  com¬ 
ponents  and  three  Reynolds  stresses  at  various  stations  in  a  180°  curved 
duct  and  downstream  tangent  of  square  cross-section.  These  are  briefly 
discussed  in  Section  4.1  of  this  report.  A  considerably  more  detailed  account 
of  this  work  is  forthcoming  and  will  be  submitted  as  part  of  a  final  report 
to  NSF  and  an  interim  (1981-1982)  report  to  ONR. 

The  quadratic  upstream  interpolation  scheme  is  also  described  in 


detail  in  an  appendix.  The  desirability  of  a  convective-differencing 
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scheme  free  of  false-dif fusion  error  is  obvious.  The  QUICK  scheme  has 
the  advantage  of  possessing  simultaneously  the  accuracy  characteristic  of 
central  differencing  and  the  stability  characteristic  of  upwind  differencing. 
As  shown  in  the  appendix  QUICK  can  be  made  to  work  efficiently  in  the 
laminar  flow  regime  and  (with  more  difficulty)  in  the  turbulent  flow  regime. 
The  higher  accuracy  per  grid  node  yielded  by  the  QUICK  approach  (relative 
to  the  standard  HYBRID  practice  embodied  in  the  TEACH  family  of  codes)*  means 
that,  in  principle,  calculations  of  a  given  accuracy  can  be  made  at  a  lower 
cost.  This  is  a  particularly  important  consideration  for  complex  3-D  flows 
(such  as  the  ones  of  interest  here)  requiring  substantial  grid  refinement 
for  their  spacial  resolution.  While  development  and  testing  of  the  QUICK 
scheme  has  been  performed  in  2-D  flow  configurations,  the  scheme  is  presently 
being  adapted  and  tested  for  the  case  of  3-D  curved  duct  flow. 


★  _  i 

From  which  the  present  numerical  procedure  ultimately  derives,  [16]. 
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4.  PRESENT  STATUS  OF  THE  RESEARCH  PROGRAM 

At  the  time  of  writing,  the  experimental,  numerical  and  turbulence 
modeling  aspects  of  this  work  continue  to  be  developed.  This  section  is 
intended  to  serve  as  a  brief  update  of  present  activities  at  the  Berkeley  campus. 
4.1  Experimental  Measurements 

Measurements  are  being  concluded  for  turbulent  flow  conditions  (Re  =  40,000) 

* 

in  a  180°  curved  duct  of  square  cross-section  corresponding  to  the  dimensions 

given  in  the  Appendix.  Two  components  of  mean  velocity  (UQ  and  Ur)  and  three 
-7  _ 

Reynolds  stresses  (u.,  u„  and  uQu  1  are  being  determined  in  detail  at 

u  r  w  r 

longitudinal  stations  corresponding  to  -5DH,  0°,  45°,  90°,  135°,  180°,  5D^, 
lOD^  and  2UDH.  Because  of  the  amount  of  work  involved,  only  the  most  critical 
longitudinal  stations  will  include  radial  component  measurements. 

Samples  of  center-plane  longitudinal  velocity  and  turbulence  intensity 
measurements  are  given  in  Figs.  1-a  and  b,  respectively.  The  data  for  0  =  45° 
was  not  available  at  the  time  of  plotting.  The  plots  are  provided  here  solely 
for  qualitative  discussion  of  the  flow.  A  report  to  be  submitted  jointly  to 
NSF  and  ONR  will  contain  a  detailed  description  of  the  flow  in  the  curved 
passage  and,  of  special  significance  to  the  ONR  program,  in  the  straight  down¬ 
stream  duct. 


This  geometry  was  originally  proposed  to  NSF  for  investigation.  The  complex 
details  of  the  flow  relaxing  in  the  straight  section  downstream  of  the  bend, 
and  the  relative  ease  and  accuracy  with  which  data  can  be  obtained  through 
flat  walls  have  made  of  this  geometry  a  critical  experimental  case  for  the 
ONR  research  program  on  turbulence  model  development. 
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The  peak  velocity  of  the  flow  at  -5DH  in  Fig.  1-a  corresponds  to  1.50 
m/s.  At  this  position  the  flow  does  not  appear  to  be  affected  by  the  bend.  At 
0°  the  center  plane  profile  has  accelerated  near  the  inner-radius  wall,  as 
would  be  expected  from  simple  potential  flow  considerations.  At  90°  and 
135°  the  main  flow  profiles  show  double  peaks  suggesting  an  intense  secondary 
motion  in  the  cross  stream  plane.  At  these  locations  the  dampening  influence 
(on  the  secondary  motion)  of  the  downstream  straight  duct  is  not  yet  felt. 

From  180°  to  20DH  the  flow  commences  a  relaxation  process  which  is  incomplete 
at  the  20Dh  downstream  station. 

Corresponding  turbulence  intensity  plots  shown  in  Fig.  1-b  point  to  various 
interesting  aspects  of  this  flow.  Like  the  mean  flow,  at  -5D^  longitudinal 
turbulence  characteristics  are  essentially  fully  developed.  At  this  station  the 
wall  values  of  turbulence  intensity  are  approximately  12%.  At  90°  and  135° 
secondary  motions  and  shearing  effects  (particularly  at  the  outer  radius  wall) 
produce. triple  peaks  in  the  turbulence  intensity  radial  distribution.  The 
implication  is  that  at  these  locations  the  flow  is  very  mixed  and  anisotropic  in 
its  characteristics.  Corresponding  evidence  of  similar  turbulence  activity  was 
not  observed  in  the  earlier  90°  bend  turbulent  flow  study  of  Humphrey,  et  al. 
[15],  It  would  seem  that  the  presence  of  the  downstream  tangent  as  of  90°  in 
[15]  had  an  earlier  stabilizing  influence  on  the  flow  in  the  bend  which,  in  the 
present  results,  appears  as  of  135°  approximately.  Similarly  to  the  mean  flow, 
the  longitudinal  intensity  profiles  show  that  the  relaxation  process  taking 
place  in  the  downstream  duct  is  far  from  having  been  completed  by  20DH. 

These  and  similar  detailed  measurements  at  other  transverse  locations 
(offset  from  the  center  plane)  will  provide  a  critical  base  case  for  testing 
the  numerical  method  and  turbulence  model  under  development. 


Figure  1-b:  Radial  distributions  of  dimensionless 
center-plane  longitudinal  turbulence  intensity, 

2  1/2 

(u0)  /Ug,  at  various  longitudinal  stations  for 
conditions  of  Figure  1-a. 
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4.2  Turbulence  Modeling 

The  details  of  the  flows  of  interest  to  this  study  are  strongly  in¬ 
fluenced  by  turbulent  transport  processes  driven  In  part,  by  large  anisotropic 
effects.  Examples  are  the  secondary  flows  arising  in  comers  driven  by 
gradients  of  turbulent  stresses.  These  can  only  be  predicted  with  a  closure 
scheme  solving  for  Reynolds  stresses  directly.  The  effective  viscosity  concept, 
based  on  the  notion  of  a  single  velocity  scale  and  length  scale,  has  been  shown 
[15]  to  yield  qualitatively  correct  predictions  of  the  principle  average 
features  of  curved  duct  flow.  However,  it  remains  to  be  explored  further. 

A  full  Reynolds  stress  closure  requires  the  solution  of  seven  transport 
equations  for  turbulence  quantities  in  a  3-D  flow.  While  obtaining  solutions 
to  such  equations  is  within  current  computational  capabilities,  there  is  a  need 
for  continued  experimentation  and  computer  optimization  to  determine  the 
empirical  constants  arising  in  these  more  complex  models.  The  simpler  approach, 
based  on  the  k/e  model  of  turbulence  requires  less  computer  time  and  core  and 
may  yield  average  flow  results  of  comparable  accuracy  to  the  higher  closure. 

The  calculations  of  [15]  were  performed  on  a  coarse  grid  and  it  is  not  certain 
at  this  time  what  refinements,  if  any,  would  improve  the  use  of  the  k/e  model 
for  predicting  the  mean  features  of  curved  duct  flows.  Both  closures  will  be 
of  interest  to  this  work,  although  emphasis  will  be  placed  on  the  Reynolds  stress 
scheme. 


These  were  fully  elliptic  calculations  and,  by  necessity,  the  grid  was 
limited  to  14  x  11  x  19  nodes.  Numerical  Impreclslons  arising  from  the 
coarse  grid  are,  therefore,  confounded  with  model  limitations.  The  semi- 
elliptic  calculation  approach  discussed  In  section  4.3  should  allow  a 
separation  of  the  two  effects. 
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In  the  two  equation  model  [17],  kinetic  energy  (k)  and  Its  dissipation  (e) 
are  determined  from  respective  transport  equations.  The  rather  simplistic 
turbulent  viscosity  concept  does  not  possess  the  required  breadth  of 
applicability  to  be  of  universal  value.  However,  turbulence  models  of  this 
kind  are  still  thriving  and  achieving  surprisingly  good  results.  One  of  the 
reasons  for  this  success  Is  well  brought  out  by  Bradshaw  [18].  In  many  flows 
the  turbulent  shear  stresses  are  at  least  two  orders  of  magnitude  smaller  than 
typical  dynamic  pressure  levels.  Therefore,  the  turbulent  shear  gradients  are 
most  significant  In  thin  localized  shear  layers.  In  flow  regions  where  stress 
gradients  are  not  large,  gradients  of  static  pressure  provide  the  main  force 
on  the  fluid  elements.  Because  turbulence  activities  are  partly  suppressed  by 
the  cross-stream  flow  at  the  walls  In  a  curved  duct,  radial  pressure  gradients 
dominate  flow  behavior.  This  point  and  economy  of  calculation  are  two  strong 
justifications  for  pursuing  a  two-equation  turbulence  model  approach. 

The  vertical  flow  emerging  from  a  curved  duct  will  commence  a  relaxation 
process  In  the  straight  downstream  section.  Every  component  of  the  Reynolds 
stress  tensor  plays  an  Important  role  in  the  relaxation  dynamics.  A  3-D 
algebraic  Reynolds  stress  model  Is  presently  being  developed  along  the  lines  of 
Launder  [19],  and  Rodl  [20]  to  model  this  complicated  flow. 

The  essence  of  the  algebraic  stress  model  lies  In  the  recognition  that 
through  various  simplifications  the  full  Reynolds  stress  equations  may  be 
reduced  to  a  set  of  algebraic  relations  among  the  Reynolds  stresses,  of  the 
general  form: 


K 


(i) 
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Appnoacnes  of  this  kind  have  been  used  by  Launder  and  Ylng  [21]  and  Tatchell 
[22]  to  predict  turbulence  driven  secondary  flows  in  ducts  of  square  cross- 
section.  Equations  of  the  form  (1)  are  obtained  if  stress  gradients  in  the 
Reynolds  stress  equations  are  eliminated.  The  convective  and  diffusive  trans¬ 
port  of  Reynolds  stresses  are  related  to  the  transport  of  turbulent  kinetic 

u1ui 

energy  by  assuming  that  the  rate  of  variation  of  — t-“-  along  a  streamline  is 
much  less  than  that  of  u.u .  Itself.  A  more  general  model  of  this  type,  in 

■  J 

principle,  should  yield  better  predictions  of  considerably  more  complex 
turbulent  flows  than  the  standard  k-e  model. 

Presently,  the  algebraic  stress  model  is  being  tested  by  reference  to 
developing  turbulent  flow  In  ducts  of  rectangular  cross-section  [23].  Sub¬ 
sequently,  predictions  of  180°  curved  duct  flow  and  of  the  flow  emerging  from 
a  180°  curved  duct  and  relaxing  in  a  straight  duct  section  will  be  compared 
with  detailed  measurements  from  this  study  (see  section  4.1)  and  two  equation 
(k/e)  model  predictions.  It  Is  anticipated  that  the  curved  and  relaxing  flow 
tests  will  be  especially  useful  in  helping  to  determine  the  algebraic  stress 
model  accuracy  and  range  of  application.  It  is  also  anticipated  that  the 
curved-straight  section  flow  comparison  will  yield  some  Insight  into  the 
manner  by  which  the  algebraic  approach  can  be  made  more  accurate  for  modeling 
complicated  3-D  anisotropic  turbulent  flows. 

4 . 3  Numerical  Procedure 

Earlier  k/e  model  calculations  of  developing  curved  duct  flows. with  strong 
curvature  effects  [15]  yielded  qualitatively  acceptable  results  but  did  not 
clarify  which  of  two  effects,  the  turbulence  model  or  numerical  diffusion, 
was  the  most  serious  limitation  in  the  calculations.  The  development  of  an 
algebraic  stress  model  of  turbulence  (section  4.2)  addresses  the  issue  of 
resolving  turbulence  model  failures.  Improvements  of  a  strictly  rwmMcal 
nature  are  discussed  here. 


In  his  doctoral  thesis,  Pratap  [13]  developed  a  semi-elliptic  calculation 
approach  which  was  applied  to  flows  in  curved  ducts  of  relatively  mild 
curvature  [6].  The  nature  of  the  scheme  is  such  that  the  pressure  field  is 
stored  three-dimensional ly  and  hence  Is  dealt  with  In  a  fully  elliptic  manner. 
However,  velocity  components  are  stored  In  two-dimensional  arrays  which, 
together  with  the  pressure  field,  are  repeatedly  iterated  upon  until  the  entire 
flow  field  satisfies  a  pre-established  convergence  criterion.  A  forward 
marching  (parabolic-type)  algorithm  Is  used  in  the  iteration  process  and 
requires  that  streamwise  diffusion  contributions  to  the  balances  of  momentum 
be  neglected.  Likewise,  regions  of  streamwise  flow  reversal  are  not  allowed. 

Notwithstanding  these  limitations,  the  semi -elliptic  scheme  is  considerably 
more  versatile  than  a  purely  parabolic  calculation  approach  and  allows  for 
economy  of  calculation  even  on  fairly  refined  grids.  It  has  the  virtue  of  not 
requiring  special  and  (generally)  complicated  approaches  for  the  calculation 
of  pressure,  such  as  those  which  are  described  in  [11,24,25], 

The  present  study  Is  concerned  mainly  with  developing  flows  of  strong 
to  mild  curvature;  2  <  RC/DH  <  20  approximately.  Parabolic  schemes  for  such 
flows  have  been  shown  to  lead  to  Inaccurate  results  [26]  and  Improvements  for 
predicting  pressure  as  outlined  in  [11,24,25]  have  only  been  tested  for  fully 
developed  flow  conditions  or  for  conditions  where  the  flow  entering  the  curved 
section  comes  directly  from  a  contraction;  the  latter  Is  a  case  of  relatively 
thin  boundary  layer  regions  and,  consequently,  localized  weak  secondary  flows. 
Fully  elliptic  calculations  [7,  15,  26]  are  expensive  and,  due  to  storage 
limitations,  cannot  usually  be  performed  with  the  degree  of  grid  refinement 
necessary  for  producing  grid- in dependent  (numerical -diffusion  uncontaminated) 
results. 
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Given  the  above,  the  semi-elliptic  calculation  approach  presents  the 
level  of  compromise  required  for  this  work  if  regions  of  possible  flow 
reversal  can  be  dealt  with  appropriately.  Small  regions  with  flow  reversal 
have  been  observed  experimentally  (and  computed)  in  the  laminar  flow  regime 
by  Humphrey,  et  al.  [7],  Measurements  available  in  the  turbulent  flow  regime 
do  not  confirm,  nor  do  they  rule  out,  the  possible  existence  of  flow  reveal, 
particularly  at  the  inner-radius  wall  of  a  curved  duct.  One  of  the  missions  of 
the  experimental  study  is  to  investigate  possible  flow  separation  and  reversal 
between  90°  and  180°  streamwise  locations. 

Calculation  instabilities  produced  by  the  (sometimes  temporary)  appearance 
of  negative  streamwise  velocities  in  a  semi-elliptic  scheme  can  be  overcome 
by  following  the  suggestion  of  Bergeles  [27]  that,  upon  calculation,  negative 
velocities  and  mass  imbalance  terms  in  the  pressure  correction  equation  be 
reset  to  zero.  Implicit  is  the  assumption  that  for  calculation  cells  within  the 
reversed  flow  region  local  continuity  is  satisfied.  While  this  practice  intro¬ 
duces  an  overall  inaccuracy  in  the  predicted  results,  the  error  is  generally 
negligible  relative  to  the  considerably  more  refined,  and  hence  more  precise, 
calculations  that  the  semi-elliptic  approach  allows.  Of  course,  in  those  cases 
where  flow  reversal  is  only  a  temporary  numerical  nuisance,  the  practice 
allows  the  calculation  scheme  to  converge  to  the  correct  final  results  without 
undesirable  numerical  instability  creeping  into  the  calculations. 

4.4  Summary 

Laminar  and  turbulent  (k/e)  calculations  using  a  semi-elliptic  procedure 
similar  to  that  described  In  [13]  have  already  been  performed  for  the  90° 
bend  flow  Investigations  of  [7]  and  [15],  Prior  to  conducting  these 
calculations  the  numerical  scheme  was  carefully  checked  by  reference  to  the 
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computation  of  developing  laminar  flow  In  a  straight  duct  of  square  cross-section 
with  and  without  one  side  wall  sliding  transverse  to  the  main  flow  direction. 
Excellent  agreement  was  found  between  the  calculations  and  the  test  case  data. 
Figure  2  shows  a  sample  of  the  results  achieved  for  the  90°  bend  laminar  flow 
case  of  [7],  Contrasted  with  the  semi-elliptic  results  are  calculations  for 
the  same  flow  using  parabolic  and  fully-elliptic  schemes.  It  is  seen  that, 
although  imperfect,  the  semi-elliptic  calculations  most  closely  approximate  the 
measurements.  Improved  agreement  was  subsequently  obtained  (using  the  semi- 
elliptic  scheme)  by  modifying  the  exit  plane  boundary  condition  treatment. 
Although  not  shown  here,  the  improved  semi-elliptic  calculation  results 
predict  very  accurately  the  peaks  and  valleys  in  the  laminar  flow  experimental 
data.  It  is  anticipated  that  implementation  of  the  higher  order  QUICK  scheme 
for  convective  differencing  will  further  enhance  the  accuracy  of  the  semi- 
elliptic  calculation  approach  for  curved  duct  flows. 

Presently,  calculations  using  the  algebraic  stress  model  are  underway 
for  the  90°  bend  turbulent  flow  of  [15]  and  the  180°  experimental  study  dis¬ 
cussed  in  section  4.1.  In  addition,  k/e  model  predictions  of  the  180°  flow  will 
be  conducted  for  comparison  with  measurements  and  algebraic  stress  results. 

In  total,  the  90°  bend  calculations  are  the  subject  for  a  communication  to 
be  presented  at  the  AFOSR-HTTfl-STANFORD  Conference  on  Complex  Turbulent  Flows: 
Comparison  of  Computation  and  Experiment.  The  purpose  of  this  conference  is 
to  evaluate  numerical  and  turbulence  modeling  schemes. 

Both  the  90°  bend  flow  calculations  and  the  180°  bend-straight  section 
measurements  and  calculations  alluded  to  in  this  report  will  be  the  subject 
of  a  final  report  to  be  submitted  jointly  to  NSF  and  ONR  in  May  1981. 
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5.  PROJECTED  RESEARCH  SCHEDULE  FOR  YEARS  TWO  AND  THREE 

This  section  provides  a  tentative  research  schedule  (in  the  form  of  a 
table)  for  work  to  be  conducted  at  Berkeley  during  years  two  and  three  of 
research.  The  schedule  reflects  the  need  to  operate  on  a  flexible  basis  with 
a  related  program  of  research  under  the  direction  of  Professor  B.E.  Launder 
at  the  University  of  Manchester  Institute  of  Science  and  Technology  (UMIST). 
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ABSTRACT 


Turbulent  flows  arise  In  numerous  engineering  applications  Involving  the 
transport  of  fluids  such  as  are  found  In  heat  exchangers,  chemical  reactors  and 
Intake  ducting  to  gas  turbines,  to  name  a  few.  Numerical  schemes  exist  which 
purport  to  model  such  flows,  and  they  require  verification.  By  using  laser 
Doppler  anemometry  (LDA) ,  experimental  data  for  turbulent  flow  In  a  sharply 
curved  duct  can  be  obtained,  which  will  help  to  verify  turbulence  models  for  this 
special  case  of  flow.  This  report  discusses  the  construction  and  Initial 
testing  of  an  experimental  apparatus  and  associated  electronic  Instrumentation 
to  be  used  In  an  on-going  Investigation  of  curved  duct  flow. 


INTRODUCTION 

Turbulent  flow  comonly  occurs  In  many  Industrial  applications.  Often  little 
detailed  knowledge  of  the  flow  Is  required,  knowledge  of  the  pressure  drop  sufficing. 
However, some  cases  such  as  Intake  ducting  to  gas  turbine  engines,  chemical  reactors 
and  heat  exchangers  require  more  detailed  Information.  Obtaining  this  Information 
experimentally  requires  considerable  time  and  expense.  Numerical  methods  for  pre¬ 
dicting  such  flows  exist  which  could  also  be  used,  at  much  less  cost.  However, 
the  reliability  of  such  methods  must  be  verlfled.and  experimental  data  which  allows 
such  verification  Is  therefore  desirable. 

This  report  describes  part  of  a  project  concerned  with  providing  such 
experimental  verification  utilizing  laser  Doppler  anenonetry  to  measure  turbulent 
flow  In  a  sharply  curved  duct  of  square  cross  section.  Specifically,  the  report 
describes  In  detail  the  experimental  apparatus,  the  LDA  system  and  the  computer 
Interface.  Preliminary  results  are  also  Included  to  show  that  the  system  as 
a  unit  functions  and  can  provide  accurate  data.  The  unit  as  a  whole  will  be 
used  In  an  on-going  Investigation  pertaining  to  turbulent  flow  In  curved  duct. 
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EXPERIMENTAL  APPARATUS 

This  section  contains  a  detailed  description  of  the  experimental  apparatus 
Including  the  flow  system,  test  section,  bulk  flow  measurement  and  valving,  and 
plumbing.  A  brief  description  of  the  system  In  operation  Is  also  Included. 

Flow  System 

The  basic  flow  system  Is  shown  In  Figures  1-a  and  b.  It  consists  of  a 
plexiglass  test  section  connected  to  an  elevated  head  tank  In  which  a  constant 
head  of  water  Is  maintained,  a  large  reservoir,  a  venturi  meter  connected  to 
a  mercury  manometer  to  measure  bulk  flow,  a  centrifugal  pump  to  feed  the  head 
tank  from  the  reservoir,  and  associated  valving  and  plumbing  trhoughout  the 
system.  Unless  otherwise  noted  the  "test  section",  shown  In  Figure  2,  consists 
of  the  whole  of  the  plexiglass  unit  Including  the  curved  section,  the  upstream 
and  downstream  tangents  and  a  flow  stralghtener  (to  be  described). 

Figures  3a  and  b  are  drawings  of  the  flow  system  with  relevant  dimensions 
Indicated.  The  reservoir  consists  of  a  4  x  4  x  4  ft  (1.2  x  1.2  x  1.2m  )  tank 
constructed  of  0.080  In.  (2.0  imi)  galvanized  steel,  coated  with  a  silicone  rubber 
sealer  for  added  protection  against  corrosion  at  the  weld  joints.  The  head  tank 
Is  2  x  2  x  2  ft  (0.6  x  0.6  x  0.6  m  )  and  of  similar  construction.  A  Unistrut 
steel  structure  supports  the  head  tank  above  the  reservoir,  with  a  distance  of 
10  ft.  4  In  (3.14  m)  from  the  bottom  of  the  reservoir  to  the  top  of  the  head 
tank. 

Pumping  Is  performed  by  a  Burkes  centrifugal  unit  (P.N.  7-G-1H.  C.N.  765-1) 
consisting  of  a  4.5  In  (11.4  cm)  Impeller  with  a  1  1/4  In  (3.18  cm)  outlet 
driven  by  a  manufacturer  fitted  3/4  HP  A.C.  motor.  The  pump  Inlet  connects 
to  the  base  of  the  reservoir  so  that  It  self-primes.  The  connection  Is  made  with 
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1.S  In  sched.  80  P.V.C.  pipe  with  two  ball  valves  between  the  pump  and  reservoir 
to  prevent  backflow  and  to  facilitate  maintenance.  Similar  piping  extends  from 
the  pump  outlet  to  a  gate  valve  operated  as  a  throttle  and  then  up  and  over  the 
lip  of  the  head  tank  and  Into  the  tank  Itself.  The  pipe  outlet  Is  horizontal. 

4  In  (10  cm)  above  the  bottom  of  the  tank,  1.5  In  (3.8  cm)  from  the  left 
wall,  4  In  (10  cm)  from  the  front  wall  and  at  a  45°  angle  to  those  walls.  This 
geometry  creates  a  gentle  three  dimensional  swirl  In  the  head  tank  allowing  the 
flow  to  rise  to  the  surface  thus  purging  any  air  entrained  before  leaving  the 
head  tank  for  the  test  section. 

Mater  discharges  from  the  head  tank  through  two  outlets  In  the  bottom. 

A  17  In  (43  cm)  section  of  2.4  In  sched.  80  P.V.C.  pipe  extends  Into  the 
head  tank  from  one  outlet  to  provide  the  constant  head.  The  other  side  of 
this  outlet  connects  to  similar  pipe  which  extends  down  to  8  In  (20  cm)  above 
the  bottom  of  the  reservoir.  A  12  In  (30.5  cm)  diameter  plexiglass  plate  held 
at  45°  from  the  bottom  of  the  tank  deflects  the  overflow  from  this  pipe  away  from 
the  pump  Inlet  to  minimize  possible  air  entrainment. 

An  8  In  (20  cm)  tall,  5.5  In  (14  cm)  Inside  diameter,  6.5  In  (17  cm)  outside 
diameter  plexiglass  cylinder, with  a  galvanized  steel  mesh  covering  the  end, 
extends  above  the  second  outlet.  This  device  serves  to  break  up  any  vortices 
which  might  entrain  air  and/or  create  a  fluctuating  head  In  the  test  section. 

A  2  In  I.D.  32 In  (81  cm)  long  brass  tube  extends  from  this  outlet  and  attaches  to 
a  2  In  I.D.  52  In  (132  cm)  long  tygon  tube  which  In  turn  connects  to  the  test 
section  entrance  flow  control  valve  (see  Figure  4). 

A  tee  ball  valve  and  a  6  In  (15  cm)  long  pipe  which  discharges  9  In  (23  cm) 
above  the  lip  of  the  reservoir  Is  plumbed  Into  the  piping  extending  from  the 
pump  to  the  head  tank.  This  valve  operated  as  a  bleed  until  a  throttle  valve 
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was  placed  on  the  pump.  It  now  serves  to  drain  the  head  tank  If  It  Is  Impossible 
to  drain  It  through  the  test  section  according  to  normal  procedure.  During  normal 
operation  this  valve  Is  closed  and  not  used. 

Test  Section 

As  stated  above,  the  test  section  consists  of  a  curved  section  with  upstream 

and  downstream  tangents  and  a  straightening  section  connected  to  the  upstream 

tangent.  The  curved  part  of  the  test  section  (see  Figure  5)  was  constructed  from 

clear  plexiglass.  One  solid  piece  3  In  (7.6  cm)  thick  was  machined  to  give  a 
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square  cross  section  channel  of  1.75  x  1.75  In  (4.45  x  4.45  cm  )  with  a  mean 
radius  of  5.875  In  (14.92  cm)  over  an  arc  of  180°.  A  plexiglass  plate  0.5  In 
(1.3  cm)  thick  was  machined  to  fit  over  the  open  channel  and  Is  held  in 
place  by  brass  machine  screws  thus  yielding  the  enclosed  curved  duct  geometry 
of  Interest  to  this  work.  An  o-rlng  seal  between  the  two  plexiglass  pieces 
ensuresthat  the  curved  test  section  Is  leak-proof  during  operation.  The 
hydraulic  diameter  DH1s  equal  to  that  of  the  straight  ducts  (1.75  In  }4.45  cm) 
and  the  curvature  ratio  Rg/0^2  *s  e9ual  *°  5*7.  The  outer  radius  wall  of 
the  curved  channel  Is  formed  by  a  0.0625  In  (1.599  mm)  thick,  curved  stainless 
steel  sheet.  There  are  four  sector  shaped  compartments  constructed  of  plexiglass 
on  the  outer  side  of  this  stainless  steel  wall.  Hot  fluid  pumped  through  these 
compartments,  or  electric  heaters.can  provide  heat  transfer  through  the  stainless 
steel  wall  to  the  flow  for  future  experiments. 

The  tangent  sections  (see  Figure  6)  are  each  55  In;  (190  cm)  long  corresponding 
to  31.4  hydraulic  diameters,  with  the  same  cross  section  as  the  curved  duct. 

With  the  flow  stralghtener  attached,  the  upstream  tangent  ensures  fully  developed 
flow  at  the  curved  duct.  The  downstream  tangent  serves  to  avoid  upstream  pro¬ 
pagation  of  disturbances  from  the  downstream  flow.  Like  the  curved  duct,  the 
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tangents  were  constructed  from  0.5  In  (1.3  cm)  thick  transparent  plexiglass. 

The  flow  straightening  section  (shown  In  Figure  5)  Is  12.5  In  (31.8  cm) 
long  corresponding  to  7.15  hydraulic  diameters.  This  section  contains  three 
removable  plates.  Each  plate  Is  a  0.125  In  (3.175  mm)  thick  piece  of  plexi¬ 
glass  containing  85  -  0.125  In  (3.175  mm)  diameter  holes  drilled  In  a  rectangular 
array  spaced  0.177  In  (4.495  mm)  on  centers  In  each  direction.  A  plate  Is  placed 
0.75  In  (1.905  cm)  from  the  Inlet  and  another  at  the  outlet  of  the  section,  with  the 
remaining  plate  In  the  middle  5.5  In  (14.0  cm)  from  each  of  the  others.  The 
perforated  plates  serve  to  homogenize  the  flow.  They  also  promote  rapid 
transition  to  a  turbulent  flow  condition  In  the  upstream  tangent.  A  5  In  (12.7  cm) 
long  1  cm  pitch  section  of  aluminum  honeycomb  placed  between  the  first  and 
second  plates  serves  to  break  up  weak  secondary  motions  Introduced  by  the  section 
of  curved  tygon  tube  between  the  constant  head  tank  and  the  test  section.  The 
entire  straightener  is  constructed  of  plexiglass. 

All  machined  plexiglass  surfaces  have  tolerances  of  ±  0.005  In  (0.127  ran). 

All  walls  of  the  tangent  sections  and  the  straightening  section  are  glued  together, 
as  are  the  heating  chamber  dividers.  The  top  walls  of  the  curved  and  straightening 
sections  are  removable,  and  fastened  with  brass  machine  screws.  0-rlngs  are  used 
to  seal  the  removable  walls  and  the  flanges  connecting  the  sections.  The  sections 
are  bolted  together  at  the  flanges  using  brass  machine  screws  and  nuts.  A  2  In 
outside  diameter,  2.5  In  (6.35  cm)  long  plexiglass  tube  Is  glued  to  the  Inlet 
side  of  the  flow  straightener  and  the  outlet  side  of  the  downstream  tangent. 

This  provides  an  attachment  point  to  connect  the  test  section  to  the  rest  of 
the  flow  system. 


Bulk  Flow  Measurement 


A  21n  Barco  venturi  meter  (P.N.  2-393),  which  Is  accurate  to  t  2%  over 
a  range  of  6-60  G.P.M.  (22.7  -  230  1/mln),  measures  the  bulk  flow.  The 
venturi  meter  Is  connected  to  a  50  In  (1.27  m)  differential  manometer  filled 
with  mercury .  The  mean  velocity  of  the  flow  In  the  test  section  is  determined 
from  the  expression: 

V  -  0.4092  [AP]1/2 

with  V  *  velocity  In  meters/second  and  AP  =  pressure  drop  In  Inches  of  mercury. 
The  Reynolds  number  may  similarly  be  determined  from  the  expression: 

Re  *  1.718  x  104  [AP]1/2  . 

Valving  and  Plumbing 

The  outlet  side  of  the  tygon  tube  from  the  head  tank  connects  to  a  brass 
adaptor  followed  by  a  2  In  bronze  gate  valve.  This  valve  serves  as  the 
throttle  for  the  test  section  and  is  followed  by  a  4  in  (10  cm)  section  of 
2  In  brass  pipe  connected  to  a  2  In  P.V.C.  union.  Another  brass  adaptor 
connects  the  union  to  a  6  In  (15  cm)  piece  of  21n  inside  diameter  tygon 
tube  attached  to  the  Inlet  plexiglass  tube  of  the  flow  stralghtener.  At  the 
outlet  side  of  the  downstream  tangent,  a  similar  tygon  tube -adaptor  union 
arrangement  Is  followed  by  a  13  in  (33  cm,6.5  D^)  section  of  2  In  diameter 
brass  pipe  connected  to  the  venturi  inlet.  The  outlet  side  of  the  venturi  is 
connected  to  a  6  In  (15  cm,  3^)  long,  2  In  brass  pipe  followed  by  a  2  In 
bronze  ball  valve  which  Is  used  when  purging  air  from  the  manometer  leads. 

An  8  In  (20  cm)  section  of  2  in  brass  pipe  follows  which  connects  to  2  In 
inside  diameter  tygon  tube  which  enters  the  reservoir  through  a  hole  In  the 
front  of  the  tank  above  the  water  level  and  then  discharges  below  the  water 
level.  This  tygon  tube  Is  securely  clamped  to  the  Unistrut  structure  supporting 


the  head  tank  to  prevent  any  oscillations  in  the  tube  which  might  propagate  back 
into  the  test  section. 


The  test  section,  plumbing,  valving  and  venturimeter  are  securely  mounted  to 
table  constructed  from  steel  unistrut.  This  table  is  securely  clamped  to  the 
floor  by  four  0.5  in.  steel  bolts  to  minimize  oscillations.  Provision  is  made 
to  level  the  table.  When  mounted  on  the  table  the  top  of  the  test  section  is 
31  In  (79  cm)  above  the  floor. 

Operation 

In  operation  mode  the  reservoir  is  filled  to  the  2  ft  (61  cm)  mark 
yielding  an  overall  head  of  7ft  9  in  (2.36  m).  Water  is  pumped  from  the 
reservoir  Into  the  head  tank  at  a  rate  of  50-60  6.P.M.  (190  -  230  tymin)  at 
wide  open  throttle.  Part  of  this  flow  (33  G.P.M.  ,  125fc/min  at  wide  open 
throttle)  empties  through  the  outlet  leading  to  the  test  section.  After 
flowing  through  the  test  section  this  water  discharges  back  into  the  reservoir 
and  continues  the  cycle.  That  portion  of  the  flow  into  the  head  tank  which 
does  not  enter  the  test  section  spills  into  the  reservoir  through  the  overflow 
outlet. 


LDA  MEASUREMENT  INSTRUMENTATION 


This  section  contains  a  description  of  the  laser-Ooppler  anemometer 
Including  mounting  and  traversing  hardware,  the  optical  components,  mirrors 
and  related  Instrumentation. 

Traversing  Mechanism 

The  laser-Doppler  anemometer  (LDA)  system  Is  shown  In  Figures  7  and  8. 

It  consists  of  a  Lexel  2  Watt  Argon  Ion  two  color  (green  -  514.5  nm  and  blue  - 
488.0  nm)  water  cooled  laser,  a  mirror  stage  (shown  In  Figures  9a,  b)  for 
reflecting  the  laser  beam  Into  the  anemometer  optics,  the  anemometer  optics 
(described  below),  and  a  4  In  (10.2  cm)  diameter  mirror  for  reflecting  the 
anemometer  beams  from  the  horizontal  to  the  vertical  direction  and  for 
collecting  backscatter  radiation  to  measure  two  velocity  components  In  the 
curved  duct. 

The  laser,  mirrors  and  optics  are  mounted  to  the  top  of  a  0.50  In  (1.27  cm) 
thick  aluminum  table  which  Is  Itself  firmly  bolted  to  an  X,Y,Z  traversing 
mechanism.  The  traversing  mechanism  is  positioned  by  three  stepping  motors 
with  linear  encoders  which  are  controlled  and  monitored  by  a  PDP  11/34 
computer.  Thus,  after  displacement  of  the  traversing  mechanism  in  any  direction 
the  new  position  of  the  measurement  volumne,  relative  to  a  fixed  reference  point. 
Is  known  to  within  5  microns.  Along  any  coordinate  direction  the  traversing 
mechanism  can  be  displaced  up  to  6  In.  (15.2  cm). 

Castors  are  mounted  on  the  bottom  of  the  traversing  mechanism  to  allow 
coarse  manual  placement  of  the  optical  Instrumentation  relative  to  the  test 
section.  Once  this  placement  Is  made,  the  unit  Is  raised  on  four  bolts  acting 
as  levelling  screw  jacks,  so  that  the  optics  are  firmly  fixed  relative  to 
the  test  section.  Fine  adjustments  with  the  stepping  motors  are  then  made  to 
obtain  the  exact  required  orientation  of  the  control  volume. 


The  LDA  optical  components,  manufactured  by  Disa  Electronics  (series  5SX) 


consist  of: 

55  x  12  beam  expander 
55  x  20  cover  and  retarder 
55  x  23  support 
55  x  25  beam  splitter 
55  x  27  beam  splitter 
55  x  29  bragg  cell* 

55  x  30  backscatter  section 
55  x  31  pinhole  section 
55  x  32  beam  translator 
55  x  34  optics 
55  x  35  color  separator 
55  x  51  300  mm  front  lens. 

Two  RCA-4526  photo  multiplier  (P.M.)  tubes  were  purchased  from  Disa  but  P.M. 
tube  housings  and  the  associated  electronics  were  constructed  and  assembled  by 
the  machine  and  electronics  shops.  One  of  these  P.M.  tube  housings  is  shown 
mounted  on  the  optics  In  Figure  7. 

The  entire  LDA  system  was  assembled  in  backscatter  mode  (the  same  optical 
components  can  be  assembled  In  forward scatter  mode).  This  arrangement  Is 
advantageous  in  systems  such  as  the  present  one  where  it  is  desirable  to  produce 
measurement  volume  displacements  without  causing  optical  misalignment.' 

The  availability  of  a  powerful  laser  source  and  the  relative  abundance  of 
adequately  scattering  particles  In  the  water  flow  facilitate  backscatter  anemometry. 
The  advantages  and  disadvantages  of  forward  and  backward  scattering  LDA  systems 
will  be  found  discussed  in  more  detail  In  [1]. 

*  In  this  work  the  Bragg  cell  frequency  shifter  was  not  used  In  conjunction  with 
the  rest  of  the  LDA  optical  components  since  reversed  flow  regions  were  not 
anticipated  at  the  location  where  measurements  were  made. 


The  characteristic  dimensions  of  the  measurement  volume  produced  by  the 
55  x  12  beam  expander  In  conjunction  with  a  310  mm  focal  length  lens,  with 
beams  Intersecting  at  a  half  angle  of  3.28°  are:  length  »  1.6  mm,  diameter 

*  100  ym.  The  number  of  fringes  produced  In  the  measurement  volume  were 
approximately  Nf  =  100.  However,  spatial  filtering  and  a  threshold  setting 
on  the  counter  reduced  the  number  of  fringes  for  signal  measurement  to 
between  10-20. 

Mirrors 

The  mirror  stage  for  reflecting  the  laser  beam  Into  the  LDA  optical 
components  Is  shown  In  Figures  7,  9a  and  9b. 

This  component  allows  a  reduction  In  the  table  length  required  to  support 
the  laser  and  LDA  optical  components  by  permitting  them  to  be  mounted  side-by- 
side  Instead  of  in  tandem.  The  mirror  stage  consists  of  three  mirror  assemblies 
mounted  to  two  parallel  17.68  x  6.5  x  0.75  in  (44.9  x  16.5  x  0.56  cm)  aluminum 
plates  held  12  In  (30.5  cm)  apart  by  2-1.5  In  (3.8  cm)  and  2-0.75  In  (1.9  cm) 
diameter  aluminum  bars.  Each  mirror  assembly  consists  of  a  Newport  Research 
Corporation  (NRC)  manufactured  1  in  (2.5  cm)  diameter  mirror*  glued  to  an 
adjustable  NRC  mirror  mount  (P.N.  MM-2)  which  is  attached  to  an  adjustable 
2.75  x  1.5  x  0.5  In  (7.0  x  3.8  x  1.3  cm)  aluminum  stand.  One  mirror  assembly 
Is  mounted  directly  to  the  lower  plate.  The  other  two  assemblies  are  bolted 
to  NRC  translation  stages  (P.N.  420-05)  which  are  both  mounted  to  the  upper  plate. 

The  horizontal  beam  from  the  laser  Is  deflected  up  90°  by  the  lower  mirror. 
The  next  mirror  deflects  this  beam  another  90°  back  to  the  horizontal,  and  the 
third  mirror  deflects  the  beam  over  another  90°  In  the  horizontal  plane. 

The  net  result  Is  a  beam  that  Is  parallel  to  the  emitted  beam  shifted  11.25  In 

*  The  mirrors  consist  of  a  pyrex  glass  substrate  with  a  dielectric  multilayer 
coating  optimized  to  reflect  the  light  at  an-lncldence  angle  of  45° 

(P.N.  10  D20  DM.2). 


(28.6  cm)  over  and  9  In  (22.9  cm)  up ,  travelling  In  the  opposite  direction. 

Various  adjustments  are  required  for  accurate  positioning  of  the  beam  along 
the  optical  axis.  A  relatively  hefty  construction  of  the  stage  unit 
Is  required  to  minimize  sensitivity  to  externally  Imposed  vibrations.  All 
aluminum  parts  were  anodized  black. 

The  anemometer  beams  leaving  the  optics  are  deflected  into  the  test  section 
by  a  4  In.  (10.2  cm)  diameter  glass-silvered  mirror  (see  Figures  7  and  8)  which 
also  reflects  the  backscatter  radiation  back  into  the  optics.  The  mirror  is 
held  In  an  NRC  post-mounted  adjustable  mirror  mount  (P.N.  625A-4),  which  Is 
supported  by  two  perpendicular  1.5  in  (3.8  cm)  diameter  aluminum  posts.  The 
vertical  post  Is  fixed  In  an  annular  receptacle  built  Into  the  optics  table, 
and  the  horizontal  bar  Is  clamped  to  it  by  means  of  an  adjustable  aluminum 
sleeve.  As  with  the  mirror  stage,  heavy  structural  components  were  utilized 
to  minimize  vibration.  All  aluminum  parts  were  anodized  black. 

Operation 

A  brief  description  of  the  operation  of  the  LDA  system  follows. 

Light  emitted  by  the  laser  is  circularly  polarized  by  the  Cover  and 
then  passes  through  the  mirror  stage.  It  then  enters  the  first  component 
of  the  optics,  the  Retarder,  where  the  circularly  polarized  light  Is  again 
linearly  polarized.  The  beam  Is  then  split  50:50  by  the  (neutral)  Beam 
Splitter  and  next  enters  the  color  Beam  Splitter  where  one  of  the  beams  Is 
separated  50:50  Into  blue  and  green  beams.  These  three  beams  then  pass 
through  the  outside  of  the  pinhole  section  and  Into  the  Translator  where 
the  relative  position  of  the  beams  can  be  adjusted  to  obtain  optimum  fringe 
spacing.  The  beams  next  pass  through  the  Expander,  allowing  the  measurement 
volume  to  be  decreased  which  Increases  the  intensity  of  the  collected  light  thereby 
improving  the  signal  to  noise  ratio.  The  beams  then  pass  through  a  convex  lens  which 


focuses  them  into  a  measurement  or  "control"  volume.  The  4  in  diameter  mirror 
deflects  these  focused  beams  Into  the  test  section. 
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The  light  scattered  by  particles  In  the  flow  Is  reflected  by  the  same  4  In 
mirror  back  through  the  lens  and  down  the  pinhole  which  acts  as  a  spatial  filter 
to  eliminate  undesirable  reflections  or  scattered  light.  This  light  Is  then 
reflected  by  a  mirror  In  the  backscatter  section  Into  the  P.M.  optics,  then 
through  a  color  separator  and  then  blue  and  green  filters.  The  separated  blue 
and  green  light  enters  respective  P.M.  tubes  where  an  electrical  signal 
Is  generated.  These  signals  are  fed  In  parallel*  to  separate  Dlsa  55196  LDA  counters 
which  validate,  filter  and  amplify  the  Doppler  signals  and  then  translate  their 
frequency  content  Into  a  digital  output  which  Is  Input  to  the  POP  11/34  (see 
next  section). 

Characteristics  of  the  Dlsa  counters  used  In  this  work  are  available  In 
[2],  For  more  detailed  descriptions  relating  to  the  use  of  counters  In  LDA 
systems  see  [3].  For  a  more  detailed  discussion  of  the  LDA  technique  see  [1]. 


In  this  work  measurements  were  only  made  of  one  component  of  velocity,  with 
the  green  (514.4  mm)  radiation. 
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DATA  ACQUISITION  AND  PROCESSING  SYSTEM 

This  section  contains  an  outline  of  the  POP  11/34  minicomputer  and  the 
data  acquisition  methodology.  For  a  more  detailed  description  of  the  computer 
system  and  data  acquisition  see  [4]. 

-  The  POP  11/34  Computer 

-  The  POP  11/34  computer,  manufactured  by  Digital  Equipment  Corporation, 
acts  as  the  central  data  acquisition  and  reduction  controller,  and  processor. 

It  has  a  34K  16  bit  word  memory  and  Is  equipped  with  dual  removable  RL01 
:d1sc  drives  (5  Mbytes  each)  and  dual  double  density  Floppy  diskette  drives 
(500  K  Mbytes  each).  The  minicomputer  Interacts  through  a  number  of  Interface 
options  with  external  devices  under  a  real  time  11  software  package  as  shown 
In  Figure  10.  Available  Input-output  devices  Include  a  Tektronix  4025  graphics 
terminal,  a  Decwrlter  II  hardcopy  terminal  and  a  Tektronix  4662  digital  pen 
plotter  (see  Figure  13). 

Data  Acquisition 

A  handshake  technique  was  employed  for  data  acquisition  whereby  a  logic 
conversion  circuit  conditioned  the  lines  from  the  01 sa  LOA  counter  to  the 
DR11-L  parallel  line  Interface  module.  At  every  validation  of  a  Doppler  burst 
a  "data  ready"  signal  was  Issued  by  the  counter  for  the  logic  conversion 
circuit.  This  circuit  then  sent  a  triggering  pulse  to  the  computer  Interface 
which  was  checked  for  data  availability  by  a  software  loop  approximately 
every  20  microseconds.  No  Interrupt  routines  were  used  for  obtaining  data 
due  to  the  higher  sample  rates  made  possible  by  the  handshake  technique. 


RESULTS  AND  DISCUSSION 


In  this  section  some  preliminary  measurements  are  provided  for  the  longitudinal 
mean  velocity  and  turbulence  Intensity  at  a  station  two  hydraulic  diameters  Into  the 
downstream  tangent  where  the  Influence  of  the  curved  section  on  the  flow  should  still 
be  noticeable.  The  main  purpose  In  making  the  measurements 
was  to  show  that  the  laser- Doppler  anemometer 

and  the  minicomputer  were  functioning  proper1y»and  to  validate  the  software  controlling 
the  entire  measurement  operation.  This  software  package  has  been  documented  In  [5] 
and  Is  briefly  outlined  below. 

The  measurement  procedure  was  at  follows.  The  optical  system  was  first  carefully 
aligned  relative  to  the  duct  measurement  station,  and  It  was  verified  that  dis¬ 
placements  of  the  LQA  by  the  traversing  mechanism  did  not  cause  misalignment.  A 
reference  point  was  then  chosen  on  one  of  the  corners  on  the  Inside  of  the  downstream 
tangent  and  was  labelled  X*0,  Y*0,  Z*0.  All  subsequent  traverslhg  was  controlled  by 
software  Instructions  relative  to  the  zeroed  point.  Instructions  were 
given  via  the  Tektronix  4025  graphics  terminal  to  the  software  package.  Indicating 
the  Initial  and  final  coordinate  positions  to  be  probed  by  the  measurement  volume, 
the  order  in  which  the  coordinate  axes  should  be  followed  and  the  AX,  AY,  AZ  1 

Increments  to  be  taken.  Various  system  constants  and  additional  parameters,  such 
as  the  number  of  data  points  to  be  taken  per  position  In  the  flow  were  also  read 
Into  the  software  package.  In  the  present  work  10  measurements  were  made  per 
position  and  each  one  of  these  was  based  on  1000  validated  Doppler  bursts. 

Before  the  experimental  run  the  counter  was  fixed  to  the  following  control 
settings:  combined  5/ 8  and  variable  fringe  number  operation  mode,  6%  comparator  accur 
low  pass  filter  set  atl  kHZ,  high  pass  filter  set  at  64  kHZ.  Optimal  conditioning 
of  the  signal  was  achieved  by  setting  the  amplifier  gain  and  threshold  window  to 
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appropriate  levels.  An  Indication  of  the  quality  of  the  Doppler  bursts  was  observed 
by  displaying  the  signals  on  a  Tektronix  5440  Scope.  Data  rateof  about  1  kHZ.wIth 
approximately  60%  val idatlon, provided  a  further  Indication  of  the  high  quality  of 
the  data. 

Upon  activation  of  the  software  package  the  minicomputer  took  control  of  the 
entire  measurement  sequence.  The  validated  data  was  temporarily  stored  In  arrays 
which  were  subsequently  processed  by  a  statistical  routine  In  the  software  package 
to  yield  a  mean  value  and  the  standard  deviation  of  the  1000  observations  per 
measurement.  Two  traverses  were  made  of  20  points  each  to  check  the  automatic 
data  aqulsltlon  sequence.  (For  a  sample  of  typical  output,  see  Table  1). 

At  each  of  the  points  In  the  traverse  10  values  of  the  mean  velocity  and 
turbulence  Intensity  were  determined.  Figures  11  and  12  show  plots  of  the 
averages  of  these  measurements.  The  size  of  the  points  provides  an  Indication 
of  the  error  In  the  measurements.  While  the  physical  meaning  of  the  results  Is 
not  entirely  clear  the  collection  of  these  results  serves  as  confirmation  of 
the  performance  of  the  automatic  data  acquisition  system.  Of  course, further 
testing  is  required  of  the  system  for  Its  complete  validation,  but  the  present 
results  are  very  encouraging. 

The  plots  shown  in  Figures  11  and  12  are  not  devoid  of  physical  meaning. 

The  profiles  were  taken  at  approximately  equal  distances  from  the  symmetry  plane 
In  the  flow.  The  precise  locations  are  shown  as  dotted  lines  In  the  squares  on 
the  respective  figures.  Inspection  of  the  nondimenslonal  velocity  measurements 
shows  the  flow  moving  more  quickly  at  high  values  of  X/D  (corresponding  to  the 
outer  radius  wall  In  the  curved  section).  The  more  striking  feature  of  the  data 
Is  contained  In  the  profiles  for  turbulence  Intensity.  These  show  high  turbulence 
activity  at  the  side  walls  of  the  downstream  tangent  but  with  higher  levels  at 


16. 


low  values  of  X/0  (corresponding  to  the  Inner  radius  wall  In  the  curved  section). 
The  latter  result  suggests  that  high  levels  of  turbulence  activity  are  being 
produced  at  the  Inner  walls  of  the  curved  section  which  are  then  convected  Into 
the  downstream  tangent  by  the  bulk  flow.  Both  figures  show  a  peak  In  the 
turbulence  Intensity  for  0.25  <  X/D  <  0.50.  Given  the  relative  precision  In 
uy  (<  5%  for  the  data  plotted)  It  Is  believed  that  the  peaks  represent  a  region 
of  strong  shear  In  the  flow.  This  Is  supported  by  corresponding  pronounced 

m 

variations  In  the  plotsfor  Uy  at  approximately  the  same  locations.  The  data  Is, 
of  course,  too  sparse  to  allow  definite  conclusions  to  be  drawn  regarding 
specific  characteristics  of  the  flow.  Further  and  considerably  more  detailed 
experimentation  Is  required. 
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CONCLUSIONS 

An  apparatus  has  been  constructed  for  the  experimental  Investigation  of 
square  cross-section  curved  duct  flow  with  long  tangents  attached. 

A  laser  Doppler  anemometer  of  Dlsa  manufacture  capable  of  making  simultaneous 
measurements  of  two  velocity  components  In  backscatter  mode  was  assembled.  The 
Instrumentation  required  various  modifications  to  make  data  aqulsltlon  feasible 
In  a  fairly  compact  space.  This  was  achieved  by  means  of  two  mirror  stages 
which  were  specifically  designed  for  this  purpose. 

The  main  features  of  the  mini-computer  -  LDA  •  automatic  data  acquisition 
system  have  been  tested  and  validated.  The  tests  Included  checking  that  the 
traversing  mechanism  responded  accurately  to  displacement  Instructions  given 
by  the  computer  and  that  data  validated  on  the  counter  was  accepted  and  sub¬ 
sequently  reduced  correctly  by  the  computer.  To  this  end  a  software  package 
described  In  [5], used  to  control  the  entire  experimental  sequence.was  debugged 

and  implemented  satisfactorily.  Further  extension  of  the  system  from  one  to  two 

* 

velocity  component  measurements  remains  to  be  dbne. 

Preliminary  measurements  at  a  station  downstream  of  the  curved  section  were 
obtained  and  are  the  basis  for  concluding  that  the  automatic  data  aqulsltlon 
system  Is  working  properly.  The  data  show  physically  plausible  and  Interesting 
characteristics.  A  detailed  Investigation  of  the  mean  and  turbulence  character¬ 
istics  of  the  flow  are  now  required. 


T 


This  Is  In  the  process  of  being  concluded. 
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Figure  la:  View  of  flow  system:  plexiglass  test  section, 
sump  and  constant  head  tanks. 


Figure  1b:  View  of  flow  system ‘.centrifugal  pump,  sump  and 
constant  head  tanks. 


Figure  2:  Transparent  experimental  test  sections:  180  degree 
curved  duct  of  square  cross-section  with  straight 
tangents  attached.  Part  of  the  laser-Doppler 
anemometer  optics  are  shown  resting  on  the  X,Y,Z 
traversing  mechanism. 
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Figure  4:  Entrance  and  exit  tangent  connections  to  the  flow 
system.  Venturi  meter  and  on-off  valve  are  shown 
connected  to  the  downstream  tangents.  Tygon  tubing 
and  a  gate-valve  are  shown  connected  to  the  flow 
straightening  section  attached  to  the  upstream 
tangent. 
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gyre  5:  Detailed  view  of  flow  test  section. 

Insert  shows  cross-section  at  90  degrees. 
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Figure  6:  Schematic  of  flow  test  section  with 
dimensions  Indicated. 
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Figure  7:  Laser- Doppler  anemometer  supported  on  motored 

and  encoded  X,Y,Z  traversing  mechanism.  Plastic 
hose  resting  on  table  provides  coolant  Mater 
to  a  2  W  Lexel  Argon  Ion  laser.  Test  section 
shown  to  left  of  photograph. 


Figure  8:  Top  end  side  view  of  laser-Doppler 
anemometer  and  X,Y,Z  traversing 
mechanism  In  relation  to  the  flow 
test  section. 
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Figure  9:  Two  views  of  180  degree  beam  deflecting 
mirror  stage. 
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Figure  10:  Laser-Doppler  veloclmeter  -  POP  11/34  Interface. 


velocity  and  turbulence  Intensity  normalized 
ilk  velocity  for  conditions  shown. 


velocity  and  turbulence  Intensity  normalized 
ilk  velocity  for  conditions  shown. 


Figure  13:  Computerized  data  acquisition  system  electronic 
Instrumentation.  From  left  to  right  are  shown: 
Tektronix  4025  Graphics  Ter.,  PDF  11/34  minicomputer, 
DISA  55L96a  counters,  frequency  shifter  domnmlxers 
and  photomultiplier  tube  power  supply,  Decwrlter 
II  hardcopy  term1na1,Tektron1x  5440  Oscilloscope. 
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Table  1:  Sample  data  output  from  automatic  data 
acquisition  system. 
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ABSTRACT 


Numerical  computations  have  been  performed  of  various  two-dimensional, 
elliptic  flows  at  high. Reynolds  number  with  a  view  to  assessing  the 
relative  merits  of  the  widely  used  hybrid  (i.e.,  upwind/ central)  interpolation 
and  the  recently  proposed  quadratic-upstream  interpolation  of  Leonard,  known 
as  QUICK.  The  latter  scheme  possesses  less  intrinsic  stability  than  hybrid 
differencing.  By  selectively  evaluating  some  of  the  convective  elements  at 
a  previous  iteration  level,  however,  satisfactorily  rapid  convergence  has 
been  obtained.  It  Is  concluded  from  the  flows  tested  that  the  extra  com¬ 
putational  effort  per  node  required  for  QUICK  over  the  hybrid  difference 
scheme  Is  more  than  repaid  by  the  greater  numerical  accuracy  that  results. 


1.  INTRODUCTION 


The  problem  of  numerical  or  “false"  diffusion  in  elliptic  flows  [1] 
plagues  the  computational  fluid  mechanicist  who  requires  accuracy  from  a 
prediction  procedure  at  a  reasonable  cost.  The  situation  is  usually  serious 
even  in  two-dimensional  flows  (where  uneconomically  fine  grids  are  often  required) 
while  it  becomes  grave  in  three-dimensional  flows  with  strong  gradients  of 
the  dependent  variables  (velocity,  pressure,  temperature,  turbulent  kinetic 
energy,  etc.).  Numerical  diffusion  also  causes  difficulties  for  the  turbulence 
modeler,  who  needs  to  eliminate  from  the  computations  any  aberrations  due  to 
the  finite-difference  practices  in  order  that  his  solutions  reflect  purely  the 
outcome  of  the  physical  model  embodied  in  the  prediction  procedure. 

In  two  recent  publications,  Leonard  [2,3]  has  reassessed  the  problem 
of  numerical  diffusion  in  convection  dominated  flows.  In  reviewing  the 
characteristics  of  the  more  commonly  used  treatments,  he  shows  [3]  how 
central  differencing  may  lead  to  unphysical  oscillatory  behavior  in  an  implicit 
scheme  or  to  non-convergence  in  an  explicit  computation  approach  in  regions 
of  the  flow  where  convection  dominates  diffusion  (Pe  >2).  He  also  shows  [3] 
that  while  upstream  differencing  of  convective  terms  leads  to  a  highly  stable 
discretization,  this  stability  is  achieved  at  the  expense  of  introducing 
truncation  error  or  numerical  diffusion  arising  from  streamline- to- grid  skew¬ 
ness,  defects  which  may  seriously  Impair  the  accuracy  of  the  computations. 

If  the  numerical  diffusion  coefficient  (r  )  is  to  be  made  insignificant  in 
comparison  with  the  physical  diffusion  coefficient,  the  mesh  Peclet  number, 
UAx/r„um  must  be  less  than  2.  While  grid  refinement  can,  in  principle, 
alleviate  the  problems  associated  with  upstream  differencing,  the  degree  of 
refinement  needed  to  eliminate  false  diffusion  is  often  impractical.  This  is 
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especially  the  case  In  problems  of  engineering  Interest,  where  although  the 
effective  turbulent  dlffuslvlty  may  be  several  orders  of  magnitude  larger  than 
the  molecular  value,  the  flow  may  be  three-dimensional,  the  boundary  con¬ 
ditions  or  flow  may  be  unsteady  and  the  field  variables  may  display  strong 
non-linear  variations.  Leonard  [2,3]  In  fact  proposed  an  alternative  dis¬ 
cretization  of  the  convective  transport  terms  which,  he  claimed,  while  greatly 
reducing  the  numerical  diffusion  created  by  an  upwind  treatment,  avoided  the 
problems  of  instabilities  found  with  central  differencing.  The  approach, 

which  Leonard,  with  a  rare  mastery  of  the  acronym  [2],  styled  tlje  QUICK 
* 

algorithm  employs  a  three-point  upstream-weighted  quadratic  interpolation. 
References  [2]  and  [3]  give  several  examples  which  demonstrate  the  Improvements 
in  accuracy  that  result  from  the  scheme.  A  very  recent  paper  by  Leschzlner 
[4]  provides  supporting  comparisons  of  the  advantages  of  both  QUICK  and 
skewed-upwind  differencing  [5]  over  the  hybrid  treatment. 

Although  these  signs  were  encouraging.  It  must  be  said  that  none  of  the 
reported  tests  have  been  numerically  very  challenging:  in  every  case  the 
equations  were  linear,  the  most  complicated  being  the  solution  of  the  thermal 
field  in  a  cavity  with  prescribed  velocity  field.  The  question  of  whether  the 
QUICK  algorithm  would  be  (or  could  be  made)  sufficiently  stable  to  justify  its 
use  in  more  realistic  cases,  where  the  flow  field  was  highly  non-linear  and 
(especially  in  the  case  of  turbulent  flow)  the  equations  intricately  coupled 
was  thus  open  to  question.  It  is  this  topic  that  the  present  contribution 


*Quadratic  upstream  interpolation  for  convective  kinematics. 
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addresses.  The  following  section  presents  the  interpolation  schemes  used  In 
this  comparative  study  while  Section  3  .discusses  the  outcome  of  tests  for 
the  test  flows  examined,  including  cases  of  both  laminar  and  turbulent  flow. 
The  paper  closes  with  some  related  conclusions  and  recommendations. 

After  the  Initial  version  of  this  paper  was  completed,  a  manuscript 
Leschzlner  and  Rodl  [6]  reporting  a  study  of  recirculating  jet  flows  came  to 
our  attention.  The  results  appearing  In  that  paper  relating  to  the  merits  of 
quadratic  differencing  are  broadly  in  accord  with  those  described  hereunder. 


2.  NUMERICAL  PROCEDURE 

2.1  General  Remarks 

The  flow  calculation  procedure  used  as  the  basic  vehicle  for  this  study 
Is  the  TEACH-2E  program  of  Gosman  and  Pun  [7].  Velocity  and  pressure  are 
solved  directly  as  primitive  variables  on  a  staggered  grid  using  the  Integrated 
control  volume  approach.  The  procedure  has  been  described  In  detail  In 
various  references;  see  for  example  Patankar  [8].  It  Incorporates  the  SIMPLE 
solution  algorithm  of  Patankar  and  Spalding  [9]  whereby,  for  a  fixed  pressure 
field,  a  set  of  difference  equations  for  the  x  and  y  velocity  components  are 
solved  line-by-line.  After  completing  a  sweep  of  the  entire  solution 
domain,  adjustments  are  made  to  the  pressure  field  to  ensure  that  continuity 
Is  satisfied  along  every  line  of  cells.  Such  adjustments  destroy  the  compliance 
of  the  velocity  and  pressure  fields  with  the  momentum  equations  and  so  further 
Iterations  are  needed  until  the  momentum  and  continuity  equations  are 
simultaneously  satisfied  to  the  required  degree  of  accuracy.  For  turbulent 
flows  the  standard  energy-dissipation  rate  eddy  viscosity  model  is  employed 
(see,  for  example.  Launder  and  Spalding  [10])  and  the  kinetic  energy  and 
dissipation  equations  are  solved  simultaneously  with  the  mean  velocity  field. 

The  two  momentum  equations  and  the  equations  expressing  the  transport 
of  turbulent  kinetic  energy,  k,  and  Its  dissipation  rate,  e,  may  be  expressed 
in  the  following  form: 

£  (Ml  (P*>  -  £  <r  $  ♦  £  <r  $  ♦  S0  0) 

where  t  Is  the  particular  variable  of  interest  (u,  v,  k,  or  e),  T  Is  the 
exchange  coefficient  and  SQ  encompasses  all  remaining  terms  and  is  therefore 
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referred  to  as  a  "source"  tern.  In  the  laminar  flow  case  (4  ■  u  and  v)  I* 

Is  slnply  the  laminar  viscosity*  y,  and  SQ  contains  only  the  pressure  gradient 
tern.  Each  equation  of  type  (1)  corresponding  to  the  different  variables  Is 
converted  to  finite-difference  form  by  Integrating  over  a  control  volume 
surrounding  each  node.  Both  the  discretization  schemes  for  convective 
transport  whose  performances  we  compare  below*  have  been  used  In  conjunction 
with  the  usual  central -difference  approximation  of  the  diffusion  terms 


(2) 


2.2  The  HYBRID  Scheme 

In  common  with  many  other  programs,  TEACH-2E  employs  basically  an  upwind 
treatment  In  representing  the  convective  terms  of  (1).  As  a  minor  Improve¬ 
ment,  however,  central  differencing  Is  used  when  the  local  cell  Reynolds 
number  Is  less  than  2.  This  hybrid  scheme  (henceforth  called  HYBRID) 
proposed  by  Spalding  [11]  may  be  stated: 
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It  Is  clear  from  equations  (4)  and  (5)  that  although  the  truncation  error  Is 
formally  proportional  to  Ax*  the  actual  size  of  the  error  depends  on  the  second 


spatial  derivative  of  $  throughout  the  computational  domain.  The  upstream 
difference  scheme  win  therefore  be  free  of  first-order  truncation  error 
(which  Is  the  main  source  of  numerical  diffusion)  only  If  the  function  t  Is 
linear  In  x;  any  function  which  has  a  local  extremum  must  Involve  a  non-zero 
first-order  truncation  error.  Leonard  [2]  argues  that,  for  convection  dominated 
flows,  HYBRID  Is  only  marginally  superior  to  pure  upwind  differencing  and  Its 
use  places  very  severe  limitations  on  the  coarseness  of  grid  that  will  allow 
useful  numerical  accuracy  to  be  obtained. 

2.3  The  QUICK  Scheme 

In  following  the  control -volume  approach  to  derive  a  set  of 
finite-difference  equations,  the  main  concern  Is  to  estimate  accurate  values 
of  the  dependent  variables  at  the  surfaces  of  the  control  volume  while  avoid¬ 
ing  the  instability  associated  with  high  convection  rates.  For  central 
differencing,  the  estimate  Is  based  on  a  linear  Interpolation  which  gives  a 
second-order  truncation  error;  but  the  scheme  becomes  unstable  for  R^  >  2. 

A  parabolic  polynomial  Interpolation,  which  simultaneously  possesses  higher 
accuracy  and  stable  convective  sensitivity,  provides  an  alternative  approach. 

The  control -volume  surface  values  are  obtained  by  fitting  a  parabola  to  the 
values  of  0  at  three  consecutive  nodal  positions,  the  two  nodes  located  on 
either  side  of  the  surface  In  question  plus  the  next  node  on  the  upstream 
side.  Examples  of  quadratic  Interpolation  expressions  for  a  uniform  grid 
spacing  are  given  by  the  following  expressions,  (see  Fig.  1): 

If  ue  >  0;  $e  »  j  +  *i+1)  -  J  (♦i_1  +  $i+1  -  2^)  (6) 

If  u0  <  0;  $e  ■  £  +  ♦•j+'j)  -  +  ^+2  "  ^1+1  ^ 


(7) 


These  expressions  may  be  Interpreted  as  linear  Interpolations  for  +e 
corrected  by  the  Inclusion  of  terms  proportional  to  the  respective  upstream- 
weighted  curvatures.  If  the  latter  terms  In  (6)  and  (7)  are  omitted,  the 
central -difference  form  Is  recovered.  Initially,  forms  (6)  and  (7)  were 
Incorporated  Into  TEACH  by  treating  the  central -difference  part  as  the  true 
convective  term  and  accommodating  the  curvature  correction  as  an  explicit 
source  term  using  the  current  value  of  the  $'s  In  question.  This  proved  to 
be  unstable  for  all  cases  examined  through  loss  of  diagonal  dominance.  False 
transient  techniques  were  evolved  that  did  succeed  In  procuring  convergence 
albeit  rather  slowly  (roughly  2%  times  as  many  Iterations  as  for  HYBRID). 
Subsequently,  several  alternative  decompositions  of  the  difference  formulae 
were  tested  and,  of  these,  the  following  alternative  groupings  secured  con¬ 
vergence  for  laminar  flows  without  the  addition  of'  false- transient  terms: 


If  ue  >  0  ;  ♦e  «  J  +  3$j+1) 

If  ue  <  0  ;  ♦g  *  ^  +  6<fri+1)  -  5  (>i+2 


source 


(8) 


If  ue  >  0  ;  *e  ■  |j  (6<>i  +  4$1+1)  -  8  ^1-1^  ^1-n) 

source 

If  ue  <  0  ;  *e  -  l  (3*1  +  4*1+1)  -  5  ($1+2  -  241+1)/ 


(9) 


source 


Equations  (8)  generally  led  to  faster  convergence  (about  35%  more  Iterations 
than  HYBRID)  but  occasionally  the  scheme  diverged;  we  would  therefore  recommend 
equations  (9)  (requiring  about  50%  more  Iterations  than  HYBRID)  as  the  safer. 

In  both  cases  the  source  parts  are  all  evaluated  with  values  determined  at 
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the  pieviouA  level  of  Iteration.  The  latter  Is  an  Important  point  to  observe 
since  we  found  that  If,  Instead,  we  used  "latest  available"  values  of  + 
convergence  was  either  destroyed  or  at  least  greatly  impaired. 

In  calculations  for  turbulent  flow,  the  strong  coupling  between  the 
turbulent  and  mean-flow  fields  has  meant  that  we  have  had  to  use  In  addition 
a  pseudo-transient  term  that  went  to  zero  as  a  converged  solution  was 
approached.  For  this  case  the  general  form  of  the  difference  equation  for 
4  at  the  center  P  of  a  cell  Is  given  by 


♦c 


E  A1  ♦i  +  $c 


Aaw,  * 

+  -sr*p 


Z  A.  +  AmP 
1  li¬ 


no) 


Convective  and  diffusive  effects,  embodied  in  the  A^  coefficients, 

are  summed  over  the  nodes  neighboring  with  P;  S  is  a  possible  source  term. 

Aiflp  * 

The  false- transient  term  (where  Am  ■  pAV  represents  the  mass  of  the 

cell  at  P,  and  *  denotes  the  value  of  4>p  at  the  previous  Iteration)  does 
not  affect  the  final  solution  although  It  significantly  affects  the  rate 
of  convergence.  No  general  rule  has  been  obtained  for  finding  the  opti¬ 
mum  value  of  At  in  the  false  transient  term  which  must  be  determined  from 
numerical  experimentation. 
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3.  TEST  CASES,  RESULTS  AND  DISCUSSION 

3.1  Axlsymmetrlc  Irrotatlonal  Stagnation  Flow 

To  obtain  some  understanding  of  the  behavior  of  the  differencing 
schemes,  the  axlsymmetrlc  Impingement  of  a  semi-infinite  extent  of  an  Ideal 
Irrotatlonal  fluid  on  a  plane  wall  was  examined.  The  enforced  turning  of  the 
streamlines  from  the  direction  normal  to  the  wall  to  that  parallel  to  It 
necessitates  that,  with  a  rectangular  grid,  the  flow  lines  cut  diagonally 
across  the  grid.  Although  not  shown  here,  the  computed  velocity  fields  for 
HYBRID*  and  QUICK  were  In  excellent  agreement  with  the  exact  solution.  How¬ 
ever,  a  significant  difference  In  the  pressure  field  was  observed  as  shown  In 
Fig.  2.  With  a  relatively  coarse  12  x  12  uniform  grid,  the  pressure 
distribution  of  QUICK  along  r/L  «  0.05  shows  very  good  agreement  with  the 
elliptic  surface  of  the  exact  solution.  The  result  for  HYBRID  differencing 
Is  far  less  accurate  with  the  same  grid  refinement.  Indeed,  even  when  the 
mesh  In  this  case  Is  refined  to  22  x  22,  differences  shown  for  the  HYBRID 
treatment  are  an  order  of  magnitude  larger  than  for  QUICK.  We  note  that 
while  the  velocity  field  Is  linear  in  this  flow  (which  Is  why  HYBRID 
differencing  leads  to  such  accurate  results  for  the  velocity  field),  the 
corresponding  pressure  distribution  Is  quadratic  In  velocity  ano  therefore 
amplifies  the  first  order  truncation  error  of  the  upstream  differencing 
approach. 


*S1nce  the  viscosity  Is  zero  In  this  problem,  the  HYBRID  treatment  Is 
Identical  to  upwind  differencing. 
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3.2  Laminar  Flow  Cases 

TWo  viscous  recirculating  flows,  the  wall-driven  cavity  and  the  flow 
downstream  from  a  sudden  pipe  expansion,  were  selected  as  a  second  class  of 
test  cases. 

Because  of  Its  geometric  simplicity  yet  strongly  elliptic  character,  the 
wall  driven  cavity  has  been  widely  used  for  examining  numerical  discretizations 
of  the  Navler-Stokes  equations.  Many  of  the  calculated  results  available 
In  the  literature  show  only  streamline  contour  plots  without  specific  velocity 
distributions.  This  has  precluded  carrying  out  direct  and  detailed  comparisons 
with  various  studies  of  relevance  to  this  work.  However,  Burggraf  [13]  has 
performed  an  excellent  numerical  study  for  the  case  of  a  square  cavity, 
using  a  sequence  of  successively  refined  grids  for  a  range  of  progressively 
Increasing  Reynolds  numbers.  He  employed  central  differences  for  the 
non-linear  convective  terms  in  the  transport  equations. 

Computations  for  HYBRID,  central  and  QUICK  schemes  were  carried  out  with 

the  TEACH  code  by  simply  replacing  the  treatment  of  non-linear  convection 

* 

terms,  respectively.  The  computed  results  of  these  three  numerical  schemes 
using  a  uniform  grid  were  compared  directly  with  velocity  distributions  obtained 
by  Burggraf  for  Reynolds  numbers  of  100  and  400  based  on  cavity  width  and  lid 
velocity.  A  comparison  among  the  schemes  Is  also  provided  at  a  Reynolds 
number  (Rfi)  of  1000. 


♦ 

A  false  transient  treatment  was  required  to  stabilize  the  central 
difference  solution. 
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At  Re  *  100,  shown  In  Fig.  3,  excellent  agreement  Is  observed  among  the 
three  schemes  for  the  variation  of  horizontal  velocity  (u)  along  the  central 
vertical  line.  At  higher  Reynolds  numbers  the  approximations  adopted  for 
convection  have  significantly  more  Influence  on  the  accuracy  of  the  final 
solution.  If  we  take  Burggraf's  results  as  correct,  at  *  400,  shown  In 
Fig.  4,  the  HYBRID  scheme  Is  less  accurate  than  the  QUICK  or  central  difference 
scheme. 

The  velocity  profile  based  on  QUICK  with  a  moderate  number  of  grid 
nodes  (22  x  22)  shows  almost  exact  agreement  with  Burggraf's  (40  x  40)  grid 
calculations.  The  present  calculations  using  central  differencing  scheme  on 
a  22  x  22  grid  show  close  but  not  as  good  agreement  with  Burggraf's  result, 
while  the  HYBRID  treatment  leads  to  a  10%  under-estimate  of  the  flow  rate  of 
recirculating  fluid. 

At  a  Reynolds  number  of  1000,  velocity  profiles  given  in  Fig.  5, 
calculated  using  the  central  and  QUICK  difference  schemes,  show  a  fairly 
close  correspondence.  Although  one  cannot  write  with  certainty,  it  would 
appear  from  the  results  at  Re  a  400  that  the  results  given  by  QUICK  would  be 
the  more  accurate.  The  HYBRID  calculations  show  a  significantly  different 
behavior,  the  recirculating  mass  flow  rate  now  being  fully  25%  less  than  for 
QUICK. 

Figure  6  compares  velocity  profiles  obtained  with  the  HYBRID  and  QUICK 
schemes  for  three  different  grid  refinements.  The  result  for  QUICK  with 
22  x  22  grid  Is  probably  accurate  enough  to  be  regarded  as  a  grid-independent 
solution  even  at  this  high  Reynolds  number.  By  contrast,  the  HYBRID  scheme 
result  for  a  30  x  30  grid  Is  fairly  similar  to  the  velocity  distribution 
generated  by  QUICK  with  an  8  x  8  grid  (which  Is  omitted  for  clarity). 

Evidently,  substantial  savings  In  computing  time  and  storage  can  be  made  with 


12 


quadratic  upstream  Interpolation,  especially  In  three-dimensional  flows. 

As  a  footnote  to  the  explorations  of  laminar  flow  In  a  driven  cavity,  we 
c 

remark  that  at  Re  ■  10  none  of  the  methods  we  have  tried  for  Incorporating 
QUICK  produced  a  converged  solution.  The  result  Is  neither  particularly 
surprising  nor  alarming  since  It  Is  unlikely  that  a  stable,  steady-flow 
solution  to  the  Navler-Stokes  equations  exist  at  such  a  high  Reynolds  number. 
In  contrast,  HYBRID  leads  unerringly  to  a  converged  result. 

The  flow  downstream  of  a  sudden  pipe  expansion  also 'provides  a  relatively 
simple  elliptic  flow.  In  contrast  to  the  driven  cavity  example,  however,  the 
shape  of  the  recirculating  region  Is  not  directly  enforced  by  the  boundary 
conditions.  For  numerical  accuracy  the  best  distribution  of  nodes  for  this 
problem  would  be  non-uniform  with  a  concentration  In  the  thin  shear  layer  that 
springs  from  the  expansion  corner.  The  present  tests  were,  however,  made  with 
a  uniform  radial  distribution  of  nodes  as  the  results  are  principally  for 
comparative  purposes.  Uniform  inlet  axial  velocity  was  prescribed,  no-slip 
conditions  on  the  walls  and  zero  streamwlse  gradient  of  u  and  v  at  the  exit. 
Figure  7  shows  axial  velocity  profiles  at  several  streamwise  positions  for  a 
Reynolds  number  based  on  inlet  velocity  and  upstream  diameter  of  500  and  a 
diametral  expansion  ratio  of  2:1.  In  contrast  to  the  cavity  flow,  there  is 
now  virtual  agreement  between  the  two  results.  The  reason  for  this  Is  that 
the  recirculation  region  is  now  highly  elongated,  reattachment  occurring 
approximately  19  step  heights  downstream.  The  streamlines  are  thus  nearly 
parallel  to  the  grid  lines  and  the  false  dlffuslvity,  which  Is  roughtly  equal 
to  0.3  [(u2  +  v2)  AxAy]**  sin  29  Is  typically  an  order  of  magnitude  less  than 
v  because  8,  the  local  angle  between  the  grid  line  and  the  velocity  vector  Is 
so  small.  Only  in  parts  of  the  recirculation  zone  and  In  the  Immediate 
vicinity  of  the  reattachment  point  do  the  streamlines  cut  at  a  substantial 
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angle  with  the  grid  and  there  velocities  are  low.  The  results  for  a  Reynolds 
number  of  1000  again  show  little  difference  between  the  two  treatments  because 
the  angle  between  the  grid  and  the  streamlines  becomes  even  smaller. 

3.3  Turbulent  Flow 

In  turbulent  flows  bounded  by  walls,  there  Is  a  very  marked  variation  of 
effective  turbulent  viscosity  over  the  flow  field.  Moreover,  the  turbulent 
viscosity  is  popularly  obtained  by  a  transport  model  in  which,  typically,  the 
local  values  of  two  scalar  properties  of  turbulence  are  obtained  from  transport 
equations  solved  simultaneously  with  those  for  the  mean  flow.  These  tur¬ 
bulent  transport  equations  contain  major  source  and  sink  terms  which  provide 
an  intricate  coupling  with  the  mean  motion  that  has  no  counterpart  In  laminar 
flow.  For  these  reasons  one  cannot  reliably  draw  conclusions  about  the 
general  efficacy  of  a  particular  discretization  scheme  from  laminar  flow  tests 
alone.  Both  the  square,  wall -driven  cavity  and  the  pipe  expansion  have  been  ex¬ 
amined  for  turbulent  flow  conditions  employing  the  standard  high  Reynolds 
number  form  of  the  k  'v*  e  viscosity  model,  [10],  For  each  case  converged 
results  were  much  harder  to  achieve  than  for  the  corresponding  laminar  flow. 

The  practice  we  have  adopted  is  to  use  the  HYBRID  solution  as  the  Initial  field 
for  beginning  the  iteration  with  QUICK.  Moreover,  the  iterations  proceed 
using  false- transient  stabilization  with  different  effective  time  steps  In 
different  regions  of  the  flow  domain.  With  these  measures,  converged 
solutions  were  typically  obtained  in  about  250  iterations.  As  with  laminar 
flow,  the  numerical  results  for  the  flow  downstream  of  a  pipe  expansion  dis¬ 
play  only  minor  differences  between  the  two  convective  approximations  and  are 
not  Included  here;  those  for  the  square  cavity  are  much  more  pronounced  and 
are  discussed  below. 


Besides  pure  HYBRID  and  QUICK  solutions,  a  nlrd  practice,  employed  by 
Leschzlner  and  Rodl  [6],  has  been  tested.  In  which  the  QUICK  approximation  is 
applied  to  the  x-  and  y-  momentum  equations  while  the  HYBRID  scheme  Is  used  for 
the  transported  turbulence  variables,  k  and  e.  The  reason  given  by  Leschzlner 
and  Rodi  for  this  mixed  treatment  was  that  the  two  turbulence  equations  were 
so  dominated  by  source  terms  that  the  values  generated  for  k  and  e  would  be 
virtually  immune  to  the  effects  of  false  diffusion.  Whether  or  not  they  had 
tried  an  orthodox  QUICK  treatment- In  the  course  of  their  work  Is  not  stated. 
Mean  velocity  distributions  along  a  central  vertical  line  (running  from  the 
floor  to  the  roof  of  the  cavity)  are  shown  in  Fig.  8  for  a  Reynolds  number  of 
10  .  Because  the  grid  Is  fairly  refined  (40  x  40)  false  diffusion  Is  small 
In  all  cases  and  differences  between  the  mean  velocity  profiles  corresponding 
to  the  different  convective  treatments  are  minor.  However,  the  coarse  grid 
(22  x  22)  calculations  of  mean  velocity  shown  in  Figure  9  display  signifi¬ 
cant  differences  among  the  various  convective  treatments.  Figure  11,  to  be 
discussed  in  detail  below,  shows  that  In  going  from  a  refined  to  a  coarse 
grid  the  level  of  turbulent  viscosity  calculated  using  the  QUICK  approach 
does  not  change  much.  The  change  is  large,  however,  for  the  HYBRID  and 
the  mixed  convective  treatments.  Thus  turbulent  viscosity  is  considerably 
underestimated  on  a  22  x  22  grid  by  the  latter  two  treatments.  It  Is 
only  the  fortuitous  compensation  by  false-diffusion  In  the  HYBRID  scheme 
which  leads  to  the  incorrect  impression  that  the  HYBRID  scheme  velocity 
profile  calculations  are  fairly  insensitive  to  changes  In  grid  refinement. 

Note  that  because  the  mixed  convective  treatment  Is  less  prone  to  numerical 
diffusion  (and  hence  Its  compensating  effects)  than  the  pure  HYBRID  approach. 
It  Is  the  mo*t  sensitive  to  changes  In  grid  refinement. 

The  distorting  effects  of  false  diffusion  are  more  vividly  displayed 
In  Figures  10-a  and  b  where  vector  velocity  plots  of  wall  driven  cavity 
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flow  corresponding  to  the  HYBRID  and  QUICK  treatments  are  contrasted. 

5 

The  calculations  were  performed  on  a  40  *  40  grid  for  a  Re  a  10  .  Clearly 

shown  for  the  QUICK  treatment  are  three  recirculation  zones,  in  the  top 

two  comers  and  the  bottom  right  hand  corner  of  the  cavity  flow  respectively, 

which  are  not  resolved  by  the  HYBRID  treatment. 

The  effect  of  grid  refinement  on  the  calculated  distribution  of 

turbulent  viscosity  is  shown  in  Fig.  11 .  Here  it  should  be  noted  that  with 

2 

the  model  adopted,  the  level  of  turbulent  viscosity  is  given  by  0.09  k  /e. 

In  the  central  region  of  the  cavity  the  levels  of  k  and  e  are  very  dependent 

upon  diffusion  since  the  generation  terms  in  their  respective  conservation  equations 

vanish  and  the  sink  terms  remain  non-zero.  The  calculated  level  of  vt  is 

thus  particularly  sensitive  to  numenLcaZ.  diffusion.  In  fact,  a  refinement  of 

the  grid  from  (22  x  22)  to  (40  x  40)  produces  roughly  a  thnzc.- hold  increase  in 

turbulent  viscosity  in  the  central  region  of  the  cavity  for  both  the  HYBRID 

and  the  mixed  treatments.  The  direction  of  the  change  is  easily  understood. 

In  the  turbulence  model,  the  diffusion  coefficient  for  turbulence  energy,  k,  is 

about  50%  greater  than  that  for  e;  each  variable,  however,  may  be  supposed 

to  be  roughly  equally  prone  to  fiaUe.  diffusion.  Consequently,  numerical 

smearing  will  lead  to  too  large  levels  of  e  and,  since  e  is  the  sink  in  the 

k  equation,  to  too  low  levels  of  turbulence  energy;  through  this  doubled-edged 

2 

effect  the  quantity  k  /e  turns  out  to  be  exceptionally  low. 

When  QUICK  is  used  for  all  dependent  variables,  the  variation  of  vt  In 
the  central  region  arising  from  the  same  grid  refinement  is  reduced  to  only 
25%.  It  Is  of  interest  to  note  that  for  the  40  x  40  grid  the  calculated  level 
of  with  QUICK  is  uniform  in  the  central  region  while  the  other  two  treatments 
display  a  waviness  similar  to  that  produced  by  QUICK  with  the  coarser  mesh. 

This  result  leaves  little  doubt  that  among  the  three  versions  examined  QUICK 
provides  the  most  accurate  representation  of  convective  transport. 


4.  CONCLUSIONS 


(I)  In  both  laminar  and  turbulent  flow  tests  Leonard's  QUICK  scheme 
has  been  found  to  give  markedly  superior  numerical  accuracy  than 
the  HYBRID  treatment,  a  conclusion  which  reinforces  and  broadens 
the  conclusions  of  Leonard  [2]  and  Leschzlner  [4]. 

(ii)  Unlike  HYBRID,  QUICK  is  not  unconditionally  stable.  Of  several 
alternative  treatments  explored  for  laminar  flow,  that  given  by 
equation  (9)  was  found  to  give  the  most  satisfactory  overall 
convergence  rates. 

(111)  For  turbulent  flow,  where  there  is  strong  non-linear  coupling 

between  the  mean-flow  and  turbulent  variables,  a  false  transient 
treatment  (equation  (10))  had  to  be  employed  to  ensure  convergence. 

(iv)  Studies  of  turbulent  flow  in  a  driven  cavity,  suggest  that,  to 
gain  the  full  benefits  of  QUICK,  the  solution  of  the  turbulent 
variables  k  and  e  as  well  as  the  mean  velocities  should  Incorpor¬ 
ate  a  quadratic  approximation  of  the  convective  terms. 
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FIGURE  CAPTIONS 


Figure  1  Nodal  configuration  for  the  east  wall  boundary  of  a  cell  at  1. 

Figure  2  Pressure  distribution  In  axlsymmetrlc,  invlscld  stagnation 
flow  along  r/L  s  0.05. 

Figure  3  Velocity  profiles  on  vertical  centreline  of  wall-driven 
square  cavity  flow  for  Re  =  100. 

Figure  4  Velocity  profiles  on  vertical  centreline  of  wall-driven  square 
cavity  flow  for  Re  =  400. 

Figure  5  Velocity  profiles  on  vertical  centreline  of  wall-driven 
square  cavity  flow  for  Re  *  1000. 

Figure  6  Effect  of  grid  refinement  on  the  calculation  of  wall-driven 
square  cavity  flow  for  Re  =  1000. 

Figure  7  Axial  velocity  profiles  for  laminar  flow  downstream  of  a 
pipe  expansion. 

Figure  8  Mean  velocity  profiles  for  turbulent  flow  In  a  square,  wall- 
driven  cavity;  calculations  on  a  40  x  40  grid. 

Figure  9  Mean  velocity  profiles  for  turbulent  flow  In  a  square,  wall- 
driven  cavity;  calculations  on  a  22  x  22  grid. 

Figure  10  Vector  velocity  plots  of  HYBRID  (A)  .and  QUICK  (B)  calculations 
of  square  wall-driven  cavity  flow  for  Re  «  105  on  a  40  x  40 
grid.  Bottom  wall  drives  the  flow,  scale  is  in  decimeters/second. 

Figure  11  Effect  of  grid  refinement  on  the  computed  distribution  of 
turbulent  viscosity. 
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